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We report on the first measurement of elliptic flow V2(pr) of multi-strange baryons S~ + S + and £2~ + Q. + in 
heavy-ion collisions. In minimum bias Au + Au collisions at % /s NN = 200 GeV, a significant amount of elliptic 
flow, comparable to other non-strange baryons, is observed for multi-strange baryons which are expected to be 
particularly sensitive to the dynamics of the partonic stage of heavy-ion collisions. The pj dependence of v% 
of the multi-strange baryons confirms the number of constituent quark scaling previously observed for lighter 
hadrons. These results support the idea that a substantial fraction of the observed collective motion is developed 
at the early partonic stage in ultra-relativistic nuclear collisions at RHIC. 

PACS numbers: 25.75.Ld 



Lattice QCD calculations, at vanishing or finite net-baryon 
density, predict a transition from the deconfined thermal- 
ized partonic matter Quark Gluon Plasma (QGP) to ordinary 
hadronic matter at a critical temperature T c w 150 — 180 MeV 
111 |2[. Measurements of hadron yields in the intermediate 
(2 < pr ^ 6 GeV/c) and high (pj > 6 — 8 GeV/c) transverse 



momentum pj region indicate that dense matter has beenpro- 
duced in Au + Au collisions at RHIC |0 H El H S 
Holl . Furthermore, previous measurements of elliptic flow 
of hadrons indicate that the matter created at RHIC is also 
strongly interacting II UM2II . Thus, in the early stage of the 
collision, dense and strongly interacting matter will lead to 
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collective effects among constituents such as transverse col- 
lective motion. If these interactions occur frequently enough, 
the system will finally reach thermalization. Due to the initial 
spatial anisotropy of the system in non-central collisions, an 
elliptic component of the collective transverse motion should 
also be present. Collectivity is cumulative throughout the 
whole collision and should survive the hadronization pro- 
cess therefore, the amount of transverse flow ob- 
served in the final state will have a contribution from the pre- 
hadronic, i.e. partonic, stage. 

Early dynamic information might be masked by later 
hadronic rescatterings. Multi-strange baryons with their l arge 
mass and presumably small hadronic cross sections HIGH 
EE El 19], should be less sensitive to hadronic rescatter- 
ing in the the later stages of the collision and therefore a 
good probe of the early stage of the collision |20]. In- 
deed, a systematic study of hadron pj spectra from high- 
energy heavy ion collisions, using a hydrodynamically in- 
spired model, shows that multi-strange baryons thermally 
freeze-out close to the point where chemical freeze-out oc- 
curs with r c h ^160 MeV Il20l 1211 which at these collision 
energies coincides with the critical temperature T c This 
may mean that multi-strange baryons are not, or much less, af- 
fected by hadronic rescatterings during the later stage of heavy 
ion collisions Their observed transverse flow would 

then primarily reflect the partonic flow. Moreover, elliptic 
flow is in itself considered to be a good tool for understanding 
the properties of the early stage of the collisions i22H23ll . pri- 
marily due to its self-quenching nature. Elliptic flow is gener- 
ated from the initial spatial anisotropy of the system created in 
non-central collisions by rescatterings among the constituents 
of the system. The generated elliptic flow will reduce the spa- 
tial anisotropy of the system and quench its own origin. Thus 
multi-strange baryon elliptic flow could be a valuable probe 
of the initial partonic system. 

In this Letter we present the first results on elliptic flow of 
multi-strange baryons E~ + and £2 + Ti^ from Au + Au 



collisions at 



: 200 GeV, as measured with the STAR de- 



tector 1 24]. About 2 million events from Au + Au collisions 
collected with a minimum bias trigger are used in this anal- 
ysis. Multi-strange baryons are reconstructed via their decay 
topology: S — > A + K and £2 — > A + K with the subsequent 
decay of A — > p + K as described in 12011 . Charged tracks 
were reconstructed in the STAR Time Projection Chamber 
(TPC) 1 25]. Simple cuts on geometry, kinematics and parti- 
cle identification via specific ionization are applied to reduce 
the combinatorial background. A detailed description of the 
analysis procedure can be found in i2olE3l . 

Figure ^ shows the invariant mass distribution for (a) 
E~ + E + and (b) £2 + £2 + candidates from minimum bias 
collisions (0-80% of the total hadronic cross-section). The 
S~ + Z + and £2 + £2 + signals appear as clear peaks around 
the rest masses (indicated by the vertical arrows) in the in- 
variant mass distribution, above a combinatorial background. 
The combinatorial background of uncorrected decay candi- 
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FIG. 1: (color online), (a) S~ + S + and (b) £l~ +Q. invariant 
mass distribution from minimum bias (0-80%) Au + Au collisions 
at y / s^"= 200 GeV. The solid line shows the combinatorial back- 
ground as estimated from a same event rotating method (see text for 
details). Azimuthal distributions with respect to the event plane of 
the (c) S~ +E and (d) Q. +£l raw yields. Dashed lines repre- 
sent the fit results. All plots shown include S~ + E + and £2~ + £2 + 
in the transverse momentum range 1 < pj < 4 GeV/c. 



dates under the peak can be determined by sampling the re- 
gions on both sides of the peak. It can also be reproduced 
by rotating the A candidates by 180° in the transverse plane 
and then reconstructing the S and £2 candidates. The rota- 
tion of the A breaks the correlation in the invariant mass and 
therefore mimics the background of uncorrected decay pairs. 
Both background determination methods provide consistent 
results. In Fig. [2(a) and (b), the combinatorial background as 
calculated from the rotation method is shown as solid lines. 
Outside the region of the corresponding mass peak, the ro- 
tation method describes the background well. The residual 
bump at lower invariant mass than the peak in Fig.Q( a ) can 
be understood as fake E candidates being reconstructed as 
Sfake(rcA,A fake (Tc random ,/? A )), where ji a and p A are the daugh- 
ters of a real A and 7t ran d m is a random Jt. The real correlation 
between Ka and p\ remains in the Efate reconstruction result- 
ing in the observed bump in the S invariant mass distribution. 
A similar mis-association happens in the £2 case with the ad- 
dition of the 7Ia being misidentified as a kaon. Our studies 
have shown that this residual correlation does not affect the 
signal peak. The raw yields are then extracted from the invari- 
ant mass distribution by counting the number of entries in the 
mass peak above the estimated background. 

The elliptic flow V2 is calculated from the distribution of 
particle raw yields as a function of azimuthal angle (j) with re- 
spect to the event plane angle v f. The S and £2 candidates 
are divided in (j) — bins, and the raw yields for each bin are 
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extracted from the invariant mass distributions as described 
above. The event plane angle *P is used as an estimate of the 
reaction plane angle i27l Eill . Here, the event plane is de- 
termined from the azimuthal distribution of charged primary 
tracks with 0.2 < pj < 2.0 GeV/c and pseudo-rapidity |r|| < 
1 .0. To avoid autocorrelations, tracks associated with a E or £2 
candidate are explicitly excluded from the event plane calcula- 
tion. Figure^shows the azimuthal distributions of raw yields 
for (c) E~ + E + and (d) £2 + £2 + with respect to the event 
plane from the minimum bias collisions in the 1 < pj < 4 
GeV/c range. To reduce the statistical uncertainties in the E 
and £2 signal extraction and because of the cos 2((j) — V P) de- 
pendence of V2, we have folded around k/2 the candidates in 
the k/2 < (j) — < K range into the n/2 > (]) — > range. 
The distributions exhibit a clear oscillation with azimuthal an- 
gle <j) — *P for both E and £2 particles indicating the presence 
of significant elliptic flow. Dashed lines are results from fit- 
ting a function gnpqn = A[l + 2v>2 cos 2((|) — V P)], where A 
is the normalization constant. Furthermore, we note that the 
amplitude of the oscillation for the E and £2 are of similar 
magnitude indicating that their V2 is similar, as will be dis- 
cussed later. The finite resolution in the event plane determi- 
nation smears out the azimuthal distributions and leads to a 
lower signal in the apparent anisotropy 1 28] . We determine 
the event plane resolution by dividing each event into random 
sub-events and determine the correction factor to be 1/0.72 
for minimum-bias collisions. In the following, all numbers 
reported on V2 are corrected for this resolution. Systematic 
uncertainties in V2 were studied by comparing the background 
determination methods described above and by changing the 
cuts used in the E and £2 reconstruction. For the E, the esti- 
mated absolute systematic uncertainties are 0.02 for the low- 
est pr bin and smaller than 0.01 for all other pj bins. For the 
£2, the absolute systematic uncertainty is 0.04 for both mea- 
sured transverse momentum bins. Correlations unrelated to 
the reaction plane (non-flow effects) can modify the appar- 
ent V2 |11]. Non-flow contributions for multi-strange baryons 
were not studied yet, but are expected to be similar to those 
calculated for A (~ -0.01 at p T =\ GeV/c and ~ -0.04 at 
p T =25 and 4.0 GeV/c) 0]. 

Figure [2] shows the results of the elliptic flow parame- 
ter V2(pt) for multi-strange baryons (a) E~ + E + and (b) 
£2 + £2 + from minimum bias (0-80%) Au + Au collisions. 
As a reference, the open symbols represent the published lHUl 
and A V2 (,£>;) from the same event class. As guideline, 
results of the fit II29I1 to V2(pr) of and A are shown as 
dot-dashed-lines. Hydrodynamic model calculations using an 
Equation Of State (EOS) with a phase transition at T { = 165 
MeV and a thermal freeze-out at T fo — 130 MeV are 
shown as dotted lines, from top to bottom, for 71, K, p, A, 
E, and £2, respectively. The expected mass ordering in hydro- 
dynamics of V2{pr) is observed with lighter particles having 
larger V2(j>t) than heavier particles. 

First, we observe in Fig. [2] (a) that for E, the V2 increases 
with pj reaching a saturation value of ~ 18% at pj ~ 3.0 
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FIG. 2: (color online). V2(pr) of (a) S _ + E and (b) £2~ + £2 
from 200 GeV Au + Au minimum bias collisions. The v>2 of and 
A 1 1 1 ] are also shown as open symbols and the results of the fits 1 29 ] 
are shown as dashed-dot lines. Hydrodynamic model calculations 
13(11 are shown as dotted curves for K, K, p, A, E~, and Q. masses, 
from top to bottom, respectively. 



GeV/c. This is similar to the result for A baryons full . In the 
lower pr region {pr < 2.5 GeV/c), the E results are in agree- 
ment with the hydrodynamic model prediction fioll . Second, 
we observe in Fig. 0(b) that the values of V2 for the £2, are 
clearly non-vanishing although they have larger statistical un- 
certainties due to their smaller abundance. Over the measured 
pj range and considering the statistical uncertainties, the V2 
of the £2 is non-zero with 99.73% confidence level (3a ef- 
fect). The £2 v'2 values are, within uncertainties, consistent 
with those measured for the E, indicating that even the triply- 
strange baryon £2 has developed significant elliptic flow in 
Au + Au collisions at RHIC. In the scenario where multi- 
strange baryons are less affected by the hadronic stage l2(ill 
and where v? develops primarily at the early stage of the col- 
lision I22H23I1 . the large V2 of multi-strange baryons reported 
in this paper shows that partonic collectivity is generated at 
RHIC. 



Previously, a particle type (baryon versus meson) difference 
in V2(pt) was observed for n and p 13111 as well as for 
and A 1 1 1] at the intermediate pj region. The present results 
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FIG. 3: (color online). Number of quark (n q ) scaled i<2 as a func- 
tion of scaled pj for S~ + S + (filled circles) and Q~ + Q. + (filled 
squares). Same distributions also shown for Jt + + ii~ (open dia- 
monds), p + p (open triangles) l3lll . (open circles), A + A (open 

squares) 1 11]. All data are from 200 GeV Au + Au minimum bias 





In summary, we reported the STAR results on multi-strange 
baryon, E~ + E + and Q. + H + , elliptic flow V2 from mini- 



mum bias Au + Au collisions at 



200 GeV. The ob- 
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collisions. The dot-dashed-line is the scaled result of the fit to K 
and A E3l. 



on the E V2(pt) follow closely the ones for A confirming that 
this observed particle type difference is a meson-baryon ef- 
fect rather than a mass effect. This particle type dependence 
of the V2(pr) is naturally accounted for by quark coalescence 
or recombination models 1321 l33l 13411 . In these hadronization 
models, hadrons are formed dominantly by coalescing mas- 
sive quarks from a partonic system with the underlying as- 
sumption of collectivity among these quarks. Should there be 
no difference in collectivity among u-, d-, and s-quarks near 
hadronization, these models predict a universal scaling of V2 
and the hadron transverse momentum pj with the number of 
constituent quarks (n q ). This scaling has previously been ob- 
served to hold within experimental uncertainties for the 
and the A when p T jn q > 0.7 GeV/c lHHl 



The n^-scaled V2 versus the n^-scaled pj are shown in Fig. [5] 
for 7T,~ + n + (open diamonds), p+p (open circles) i3lll . Kg 
(open triangles), A + A(open squares) fiH, E~+E + (filled 
circles) and £2 + Q. + (filled squares). Except for pions, all 
hadrons including E and £2 scale well within statistics. The 
discrepancy in the pion vi may in part be attributed to its 
Goldstone boson nature (its mass is smaller than the sum of 
its constituent quark masses) or to the effects of resonance de- 
cays (a large fraction of the measured pions will come from 
the decays of resonances at higher pj) 11291 13511 . This further 
success of the coalescence models in describing the multi- 
strange baryon \'2{pr) also lends strong support to the find- 
ing that collectivity developed in the partonic stage at RHIC. 
In addition, the good agreement of \>2(pt / n q ) /n q for p(uud), 
A(uds), Z(dss) and Cl(sss) further supports the idea that the 
partonic flow of s quarks is similar to that of u,d quarks. Fu- 
ture measurements with higher statistics, specially for the £2, 
will allow for a more quantitative comparison. 



servations of sizable elliptic flow and the constituent quark 
scaling behavior for the multi-strange baryons suggest that 
substantial collective motion has been developed prior to 
hadronization in the high-energy nuclear collisions at RHIC. 
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